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Abstract 

After setting up a Hamiltonian formulation of planar (matrix) quantum mechanics, we 
illustrate its effectiveness in a non-trivial supersymmetric example. The numerical and 
analytical study of two sectors of the model, as a function of ’t Hooft’s coupling A, reveals 
both a phase transition at A = 1 (disappearence of the mass gap and discontinuous jump 
in Witten’s index) and a new form of strong-weak duality for A —► 1/A. 
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1 Introduction 


More than thirty years after its introduction jTj the large -N limit of four-dimensional QCD 
remains elusive. It is widely believed that such a limit should capture the most interesting 
non-perturbative properties of QCD, such as confinement and spontaneous chiral symetry 
breaking, while neglecting others (e.g. resonance widths and the U( 1) anomaly). It is 
also believed that, in’t Hooft’s limit, QCD should lend itself to an effective description in 
terms of a string theory (see e.g. [2] for recent developments of the subject) or, perhaps, 
of a gravitational dual similar to the one enjoyed by J\f = 4 super-Yang-Mills (SYM) 
theory through the AdS-CFT correspondence [3]. Indeed, the large-iV classification of 
diagrams according to topology closely resembles the loop expansion of string theory in 
terms of surfaces of increasing genus. 

One should keep in mind, however, that the connection between large-IV and graph- 
topology is only proven order-by-order in perturbation theory. Whether the true large -N 
limit (defined as solving exactly the theory at finite N and then taking N to infinity) 
does actually coincide with the non-perturbative solution of a suitably defined “lowest 
genus” theory remains to be proven case by case. Already the Gross-Witten model [3j 
and the work of Marinari and Parisi |5J have taught us that the N —> oo limit may not 
commute with other limits, such as the full resummation of the strong-coupling expansion 
or approaching first a phase transition. In QCD itself, the assumption (now supported by 
lattice calculations 0) that the topological susceptibility depends on whether the large -N 
limit is taken before or after the chiral limit provides a solution of the 1/(1) problem {7j. 

In this letter we shall consider a class of planar matrix models in a (possibly new) 
Hamiltonian formulation. (Large -N literature being very vast, we shall refer to the nice 
reprint collection |3j| for the classic papers on the subject). We shall then solve a particular 
supersymmetric case (both numerically and analytically) in a planar (or better lowest- 
genus) approximation and point out several amusing features of the solution, including a 
phase transition and a non-trivial strong-weak duality in’t Hooft’s coupling A. 

One motivation for this work was to prepare the ground for checking, in a simpler 
context, a recently claimed planar equivalence JSj between a supersymmetric “parent 
theory” and its non-supersymmetric “daughter” in a particular subsector. Our method 
(or at least its numerical part) should apply without major modifications to the latter 
theory, and therefore such a check should be forthcoming. Another motivation came 
from the recent studies of the supersymmetric Yang-Mills quantum mechanics in various 
dimensions ma Although done at present mainly for SU(2) [TT] (and partly for SU(3) 
and SU(4) [T2J) gauge groups, the goal of these works is to extrapolate eventually towards 
large N whereby making contact with M-theory na. Our results here offer the prospect 
of a substantial shortcut for the whole program. 

The rest of this letter (see !2| for more detailed account) is organized as follows: we 
first define a general class of planar quantum mechanics (PQM) models and specify some 
conditions for our method to be applicable. We then focus on the supersymmetric case 
and, eventually, on a particular example for which spectrum and main features can be 
worked out both numerically and analytically. 
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2 Hamiltonian formulation of PQM 

Our formulation of planar quantum mechanics (PQM) is best done directly in a standard 
Hamiltonian framework (see also ra ). Let us start by defining a Hilbert space (with 
states that span it) and operators (acting on it). The operators will be N x N destruction 
and creation matrices: 

= ; k = l,...N f , (1) 

where i,j represent “colour” indices while k represents a generalized “flavour” index. The 
latter can be used, in a QFT generalization, as a label for momenta, polarizations etc. In 
this paper it will only serve the purpose of distinguishing bosonic and fermionic degrees 
of freedom. The above operators are assumed to obey standard (anti) commutation 
relations. In familiar notations: 

WM^}=5™5i5 k , k ,. ( 2 ) 

The Hilbert space is constructed out of the usual Fock vacuum (annihilated by all 
( fc )j\T/) by acting on it with a single-trace string of creation operators. This is the first 
essential difference between general and planar QM. It is of course mimiking the colour 
structure of the states that propagate in genus-zero diagrams, a small subset of all states 
that are singlets of a U(N) group acting as: 

M -»• UMU ] . (3) 

This is of course an enormous simplification: if, for instance, the index k takes a single 
(bosonic) value, the states spanning the Hilbert space are just labelled by a single integer, 
n, corresponding to the number of Mbs in the trace (see below). 

The Hamiltonians we shall consider are also single-trace operators and therefore sin¬ 
glets of U(N ), but, of course, they will include both creation and destruction operators. 
In order to avoid producing unwanted vacuum diagrams we shall impose that the Hamil¬ 
tonian H annihilates the Fock vacuum. This typically implies (though it is not implied 
by) normal ordering of the operators appearing in H (not to be confused, of course, with 
the order in the colour trace). 

The Hamiltonian will be a sum of such single-trace operators and will contain a factor 
g n ~ 2 for a term containing n operators. Schematically: 

H = Y J Cng n - 2 Tr(M n ), (4) 

n 

where M n stands for a product of Ms and M^s with a total of n factors. The (’t Hooft) 
limit to be consider is, as usual, N —> oo with the’t Hooft coupling A = g 2 N kept fixed. 

When such a Hamiltonian acts on a generic single-trace state it will not give, gener- 
ically, another single-trace state. However, whenever it does not, one gets subleading 
terms in the large- N limit. If we discard such terms we have a closed system and the 
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matrix elements of the Hamiltonian turn out to be functions of A alone: we simply have 
a well-defined Hamiltonian to diagonalize in the single-trace Hilbert space. 

The final ingredient of our approach is to introduce a cut-off B = n max in the occu¬ 
pation number thus reducing the problem to one that can be managed numerically P! . 
Eventually, by increasing B , we can check whether the lowest eigenstates and eigenvalues 
converge to some finite limit. As we shall see, this will be the case in a simple toy model 
where the method gives very interesting indications of the dependence of the spectrum 
from A. In turn, the numerical results will suggest properties that we shall be able to 
derive analytically. 

3 A class of supersymmetric matrix models 

We will now specialize to the case in which there is just one bosonic and one fermionic 
matrix, denoted, respectively, by a and / (plus their Hermitian conjugates). 

The class of supersymmetric models that we consider are a straightforward matrix 
generalization |5j of Witten’s supersymmetric quantum mechanics (SQM)JTB}. We will 
assume that the supersymmetric charges Q and Q f are single-trace operators that are 
linear in the fermionic matrices / and /h Thus: 

Q = Tr(A t f) , Q f = Tr(Ap) , (5) 

where A = A(a, cd) represents some function of the bosonic matrices. We also demand 
that Q and Q t are nihlpotent. This gives: 

[A®, Af\ = 0 , for all i,j , A;, l. (6) 

In our explicit example we will satisfy this condition trivially by making A (kb) depend 
only on a (cd). We then obtain the supersymmetric Hamiltonian as: 

H = {Q',Q} = (/*)$ ft[A{, A) + V' . (7) 

By construction, H commutes with Q and Q'. It also commutes with fermion number, 
F = Tr(/t /). In order to get rid of disconnected diagrams we need the condition: 

A* 10 >= 0 , for all i , j , (8) 

which is again trivially satisfied if A depends only on a, something that guarantees that 
also Q and Q ! annihilate the trivial (empty) Fock state. Thus, by construction, our model 
has (at least) one zero-energy state and does not break supersymmetry. 

The spectrum of the theory should then consist of a zero-energy sector (providing a 
certain value of Witten’s index) and degenerate massive supermultiplets. It is easy to 
show 1 that, barring unexpected extra symmetries, these supermultiplets should consist of 
just one boson and one fermion. Technically, this is a consequence of the fact that the 

1 One of us (GV) wishes to thank A. Veinshtein for a useful discussion on this issue. 
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only other operator in the SUSY algebra (besides Q, and H), C = [Q\Q], satisfies 
the equation C 2 = H 2 . Hence, eigenstates can be classified according to whether they 
carry C = ±H = ±E (in amusing analogy with BPS states). Furthermore, the algebra 
implies that states with positive (negative) C are annihilated by Q * (Q), while they are 
transformed in a state with opposite C by the other supersymmetric charge. All non¬ 
zero-energy levels must therefore consist of two states with opposite C'-parity. For the 
F — 0 (F — 1) sector C is negative (positive) for all the states but this fails to be the 
case for higher values of F (see Section 6). 

4 A specific model and its numerical analysis 

We now specialize further our model by taking: 

Q = Tr[fa\ 1 + get)] = Tr[fA\ Qt = Tr[/t(i + ga)a] = Tr[pA], (9) 

and therefore 

h = h b + h f - (10) 

H b = Tr[a^a + g(a^a + a)a 2 ) + g 2 a^a 2 ] ; (11) 

H F = Tr[/ t / + g(Pf(a f + a) + / f (a f + a)/) 

+ fi ,2 (/ t a/a t + /W/ + / t /a t a + f ] a ] fa)}. (12) 

In most of this paper we shall limit our attention to the F = 0 and F — 1 sectors. Some 
discussion of our expectations for the F > 2 will be given at the end of the paper but a 
detailed analysis is postponed to further work. 

As already anticipated, the planar states in the F = 0 sector are simply labeled by 
the integer n corresponding to the number of a's in the trace. We shall denote the 
normalized state with n bosonic quanta by |0,n). Similarly, the generic (single-trace) 
F — 1 normalized state will be denoted by |l,n — 1): the corresponding creation operator 
contains one fermionic and n — 1 bosonic operators. In the free theory (g = 0) there is a 
single zero-energy bosonic state, |0,0), while |0,n) and |1 ,n— 1) form a supermultiplet. 

Working out the matrix elements of the Hamiltonian is straightforward although te¬ 
dious (in particular the normalization factors have to be kept accurately). Keeping only 
the leading terms as N —> oo we find that the final result for the matrix elements of H 


depend only on A = g 2 N, and are given by: 

(0, n\H\0, n) = (1 + A(1 - 5 nl )) n , (13) 

(0,n + l\H\0,n) = (0,n\H\0,n+ 1) = V\n(n + 1), (14) 

(l,n\H\l,n) = (■ n + 1)(1 + A) + A , (15) 

(1, n + 11iY| 1, n) = (1, n\H\l, n + 1) = \/A(2 + n). (16) 
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After introducing a cutoff B = n max we can diagonalize the Hamiltonians in the two 
sectors and compute the spectra. Eigenvalues with E « B converge rapidly to some 
finite values, except near A = 1, where some critical slowdown of the calculation emerges. 
Fig. 1 gives the lowest fermionic and bosonic eigenvalues as functions of A. Apart from 
the trivial bosonic ground state we observe that: 

• There is excellent boson-fermion degeneracy if we stay away from A = 1; 

• The cutoff B explicitly breaks supersymmetry, which we expect to recover only at 
B = oo; 

• The breaking of SUSY allows the supermultiplets to split near A = 1. More amus¬ 
ingly, above A = 1, the supermultiplets form once more, but with new partners. The 
10, Ei) state remains unpaired (with zero energy), while \0,E n+2 ) pairs with |1 ,E n ) 
(n — 0,1,,. ..) rather than with the small-coupling partner 11, E n+ 1 ); 

• Eigenvalues tend to collapse to zero at A = 1 as the cutoff is increased. 

• Some kind of symmetry between strong and weak ’t Hooft coupling appears. 

Obviously, the behaviour near A = 1 is strongly suggestive of a phase transition (if 
the cutoff is removed). A rather shocking way of showing this is to plot the Witten index 
and partition function (restricted to the two sectors we have considered): 

W(l3,\) = Tr((-l) F e- fH ) , Z((3, A) = Tr (e~ ,3H ) (17) 

The results are shown in Fig. 2. The sudden jump by one unit in W(f3, A) around A = 1 
is quite spectacular. The standard, thermal partition function shows even more dramatic 
singularity at A = 1. As expected, the large cutoff and large (3 limits do not commute. 
Our numerical results suggest that the large cutoff limit at fixed (3 reveals a singularity of 
Z at A = 1 which tends to a step-function (i.e. as for W) plus a kind of “delta function” 
as (3 —> oo. 

In order to understand better these numerical results and what they mean at infinite 
cutoff, we now resort to some analytic methods. 

5 Analytic solution 

Let us introduce new “composite” creation and annihilation operators for single trace 
states: 

a' n ( a n ) creates (annihilates) |0, n) ; , (/ n ) creates (annihilates) |l,n-l), (18) 

that (anti)commute as usual. Introducing for convenience b = V% it is easy to show that: 

OO OO 

//U=°) _ a t l0l + y Al + 6 2 )a t n a n + (y b\/n(n + 1 )a l n a n+1 + h.c .) (19) 

n =2 n =1 


6 







2 


40 £ B £ 200 


10 £ B £ 200 


WH6,1L 




Figure 2: A dependence of the Witten index and the partition functiom, at /3 = 6 for 
different cutoffs. 

and 

OO OO 

H {F=1) =Yl[ n+ ( n+ 1 ) fe2 ]/n/n + 1 + h.C.) (20) 

n— 1 n— 1 

We can also construct the SUSY charges as: 

OO 

Q — a \fi + \/n + la^ n+ i(/ n+ i + bf n ), (21) 

n=l 

(and similarly for Ql) and check that {Q, Qt} = H. 

In the F — 0 sector, besides the trivial vacuum, |0)i, one can formally construct a 
second state annihilated by H: 

CO / \ n j 

|»>^ T < 22 ) 

Clearly its norm is only hnite at b > 1 explaining why there is no such a zero-energy state 
below b = 1. Using the formula given above for one can also check that Q' 10 )2 = 0. 

We now come to the derivation of a new sort of strong-weak-A duality which surpris¬ 
ingly exists in this model. Using (EH for both b and 1/6 we find immediately: 

bH {F=1 \l/b) = -H [F=1 \b ) + (1/6 - 6), (23) 

6 

sillCe En fnfn = 1 in this sector. Because of SUSY it must also work in the F = 0 sector. 
In terms of eigenenergies, the duality relations read: 

b (Bf-‘)(l/6) + 1) = 1 (£f =1) (6) + 1) ; b (Bf-°>(l/6) + l) = \ (b£t 0) (6) + l) . 


(24) 








b 

1/2 

2 

3/4 

4/3 

n B 

10 

3.18818 

3.19736 

1.68608 

1.91640 

1 20 

3.18808 

3.18808 

1.67497 

1.68238 

30 

3.18808 

3.18808 

1.67488 

1.67499 

exact 

3.18807663 

1.67488116 

20 

6.45527 

6.45637 

3.09853 

3.45017 

3 40 

6.45524 

6.45524 

2.97193 

2.97543 

60 

6.45524 

6.45524 

2.97177 

2.97177 

exact 

6.45523985 

2.97177086 

30 

9.49504 

9.49513 

4.46674 

4.84625 

5 50 

9.49503 

9.49503 

4.19908 

4.21714 

70 

9.49503 

9.49503 

4.19672 

4.19675 

exact 

9.49503451 

4.19671404 


Table 1: The cutoff dependence of three eigenenergies at two pairs of b values related by 
the duality described in the text. Exact results (c.f. below) are identical for dual partners. 

Notice that, due to the existence of the second vacuum for b > 1, states, in the F — 0 
sector, whose energies are related by duality, do not have the same n. These duality 
relations are nicely satisfied by our numerical eigenvalues computed at large cutoff (see 
Table 1). 

A consequence of duality is that, for all levels, the left and right derivatives of E at 
A = 1 (E'^E'y) should satisfy : 

E , > + E , < = 1/2 (25) 

a relation that has also been checked numerically. 

We finally turn to an analytic determination of the massive spectrum. To this purpose 
it is convenient to rewrite El in the F — 0 subspace as: 

OO 

H = ^2 B i B n , Bn = Vna n + bVn + la n+ i, 

n— 1 


and to introduce new states: 

| B n ) = Bl |0) = \fn\n) + b\Jn + 1| n + 1). (26) 

These states are not orthonormal, nevertheless they form a complete set and this suffices 
for our construction. The action of a “reduced” Hamiltonian, H = (El — b 2 )/b, on the 
| B n > states is very simple: 

H\B n ) = n\B n _i) + n (b + —^ | B n ) + ( n + l)|H ra+ i), n — 2,.... (27) 
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with two exceptions at n = 0,1 for which: 

H\B 0 ) = -b\B 0 ) + IB,), H\B i> = (b + ^ | B,) + 2|H 2 ). 

The simplicity of (EZD allows us to map the eigenproblem of H into a simple differential 
equation. Let us expand any generic eigenstate of H into the | B n ) basis and associate 
with it a function of one variable x: 

OO OO 

\i>) = ^2cn\B n ) , f(x) = ^c n x n , (28) 

n= 0 n=0 

which is in fact a generating function for the {c n } coefficients. 

It is then easy to see that the eigenequation for \i[)) maps into the following first order 
differential equation for f(x): 

w(x)f'(x) + xf(x) - 6/(0) - /'(0) = ef (x), (29) 


where w(x) = (x + b)(x + 1/b), and e is the eigenvalue of H (E — b(e + b)). The solution 
of (12911 is straightforward: 


where: 


f(x) = g(x) 



b_m±m_ , 

w(x')g(x') 


g(x ) — (x + b) a (x + 1/6) 


a— 1 


a = 


e + 6 
6 - 1/6 ’ 


E = a(b 2 - 1), 


(30) 


(31) 


is a solution of the homogenous equation and Xq is to be determined by some boundary 
condition. 

However, since the inhomogeneous term is given by 6/(0) + /'(0), there is an additional 
consistency condition, namely the solution and its derivative, when taken at x — 0, must 
reproduce /(0) and /'(0) again. This leads to the relations 


either (6p(0) + </(0)) = 0, or / dx(x + b) a 1 (x+l/b) " = 0. (32) 

J Xo 

The first condition gives a = 0, hence it can only lead to zero-energy states. Thus the 
massive spectrum follows from the second condition 432D- Consistency with (1221) fixes 
x 0 = —1/6 for 6 > 1 and x 0 = —6 for 6 < 1. In either case one should use the analytic 
continuations of Eq. (EHD in order to solve (E2D- 

Once this is done, our solution can be written in terms of the standard hypergeometric 
function F = 2 F, as 

/w = S^TI76 F(1 ' Q;1 + a; ^n76 ) ' 6<1 ’ (33) 

/W = -l- 0(U-a;2-«ii±/), b> 1, (34) 

1 — a x + b x + b 
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and provides the generating functions for the expansion coefficients {c n } of the arbitrary 
eigenstate into the | B n ) basis. As one crosscheck examine the b > 1 solution at a = 0 
to find that it indeed generates the second massless state, Ea. (l221l . On the other hand, 
similar attempt for the b < 1 solution fails - there is no such state in the weak coupling 
regime. 

To summarize, after some trivial change of integration variable, the non-zero-energy 
levels of the F = 0 (and thus by SUSY also of the F — 1) sector are given by the roots 
in a of the following equations: 

/•1/A r\ 

/ dx(l — = 0 , (A > 1) ; / dx(l — x)~ 1 x a ~ 1 = 0 , (A < 1). (35) 

Jo Jo 

Solving these equations indeed reveals a series of discrete zeros, a n > (<)0 for b > (<)1 
which nicely confirm the eigenvalues E n = a n (A —1) computed numerically in the previous 
section. One immediately checks that, for two values of A related by A —> 1/A, the 
solutions for a are connected by a —> 1 — cc, insuring the duality relations fl2H) among the 
corresponding eigenvalues. 

At this point one can study the flow of the eigenvalues in various situations, e.g. at very 
weak (and thus by duality also at very strong) coupling. More interesting is the behaviour 
of the eigenvalues near the critical point at A = 1. The Beta-functions in (E3) can be 
related to F(a, 1, a + 1; A) and F( 1 — a, 1,2 — a; 1/A) for A < 1 and A > 1, respectively. 
^Frorn the known asymptotic behaviour of F (as its last argument approaches 1), we 
easily get the approximate eigenvalues around A = 1 in the form: 

A -> 1“ : E n = (-a n )( 1 — A) , 

A -> 1 + : E n+ 1 = (1 - a n ) (A - 1) , 7i = 0,1, 2,... , 

fp(oin) + 7 + l °g(\l - A|) + O (|1 - X\log(\l - A|)) = 0 , (36) 

where is the logarithmic derivative of the T-function and 7 = 0.5772.. is the Euler- 
Mascheroni constant. These formulas obey the duality relations (P-lh They also show 
the non-analytic way the various levels collapse to zero energy at the critical point. In 
particular, as A —> 1 _ , the first eigenvalue approaches zero as —(1 — A) log~ l (l — A), i.e. 
with vanishing first -and infinite second- derivative. 

The above formulae also allow a quantitative study of the free energy of the model 
near the phase transition, which appears to be stronger than in the Gross-Witten model. 
We have also determined numerically the first few zeroes at generic values of A, and found 
perfect agreement with the large cutoff limits of the numerical eigenvalues. The only 
slightly difficult comparison occurs just around the phase transition where convergence 
(as one increases the cutoff) undergoes a critical slowdown. 

6 Remarks about the F > 2 sectors 

In principle our analysis can be extended in a straightforward way to higher fermion- 
number sectors. In practice, calculation of the Hamiltonian in those sectors becomes 
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Figure 3: Chew-Frautschi-like plot of the lowest seven levels with their respective degen¬ 
eracy and supersymmetry partnership at vanishing A. 
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quickly cumbersome. There are at least two reasons why it would be worthwhile making 
such an extension. 

Firstly, one would like to check whether the F — 0,1 sectors completely determine the 
structure of the phase transition at A = 1. This would depend on how eigenvalues in the 
higher -F sectors behave near A = 1 and in particular on whether there are discontinuous 
jumps of Witten’s index also in those sectors. 

The second reason is that we may expect qualitatively new phenomena to occur when 
we consider F > 2 sectors. The fact that eigenstates in the F = 0 and F — 1 sectors pair 
nicely without involving, say, F = 2 states can be argued on the basis of simply counting 
the former states at weak coupling. However, when we go to higher F, states are typically 
highly degenerate at zero coupling. The counting is relatively easy and is summarized in 
a kind of “Chew-Frautschi” plot in Fig. 3, where non-degenerate states are marked with 
a full circle while degenerate ones are represented by an open circle showing the degree 
of degeneracy. 

Pairing these states in SUSY doublets (as A is switched on) turns out to be possible, 
but non trivial (due to some magic combinatorics! jQHj-IEI): it is shown in Fig. 3 via the 
vertical segments connecting different circles. For instance, the E = 6 levels must pair 
according to the following pattern: the two F = 2 states find their SUSY partners in two 
linear combinations of the four F — 3 states. The remaining two F = 3 states will match 
two (linear combinations) of the three F — 4 bosons. Finally, the third F — 4 boson will 
pair with the single F = 5 fermion. One can also argue that, while most of the above 
levels have C/E = (—1) F+1 , two of the F = 3, E = 6 levels and one of the F = 4, E = 6 
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levels have C/E = (—1) F . We are planning to check numerically the low-lying spectra of 
the F — 2 and F — 3 sectors in order to see whether, indeed, two of the four F — 3 states 
around E — 6 do split from the two (F — 2, F — 3) doublets as we turn on A. Were this 
not the case, would signal some higher symmetry underlying the model. 

7 Discussion, summary 

In this paper we have presented a new way to tackle, both numerically and analytically, 
planar quantum mechanical problems which hopefully represent the large -N limit of ma¬ 
trix models. Given the ubiquiness of matrix models in theoretical physics (see again 0 ), it 
is hard to overestimate the importance of developing powerful techniques for approaching 
this kind of questions. 

Our method is based on a direct Hamiltonian construction of states and operators that 
are relevant at lowest genus in a topological expansion of the theory. In principle it should 
be applicable to any discretization of quantum field theories that allows to compute the 
planar Hamiltonian in a convenient basis of vectors. 

As an illustration of the method we have considered a supersymmetric quantum me¬ 
chanics model and managed to solve for its spectrum, both numerically and analytically, 
in two fermionic sectors. Since supersymmetry transformations close within these two 
sectors, we find, as expected, boson-fermion degeneracy. To our surprise, however, we 
also find that, at a critical value of the ’t Hooft coupling, A = 1, the spectrum looses 
its mass gap and becomes continuous. This conclusion is also confirmed by the radically 
different dependence of the spectrum on the cutoff at A = 1. This dependence is indeed 
characteristic of the scattering plane-waves El. On the other side of the critical value 
the spectrum has once more a mass gap but there is one more zero-energy bosonic state. 
In other words the Witten index has jumped by one unit across the phase transition. 
Furthermore, energy levels on the two sides of the critical point are connected through 
a non-trivial duality relation. Another amusing property of the model is that, at least 
within those two sectors, it can be solved analytically. 

Besides generalizing the model to more interesting cases, there are two important 
directions in which the model itself deserves further study: 

• Extend calculations to sectors with higher fermion number; 

• Understand the situation at finite (though large) N. 

We hope to be able to address these issues in a forthcoming paper. 
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